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Abstract 



We derive contributions to the trace formula for the spectral density accounting for the role of 
^ ■ diffractive orbits in two-dimensional billiard systems with corners. This is achieved by using the 
', exact Sommerfeld solution for the Green function of a wedge. We obtain a uniformly valid formula 
" which interpolates between formerly separate approaches (the geometrical theory of diffraction and 
Gutzwiller's trace formula). It yields excellent numerical agreement with exact quantum results, also 
_ . in cases where other methods fail. 
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1 Introduction 



Two-dimensional classical billiards became popular as model systems exhibiting a rich variety of 
dynamical behaviour, ranging from integrable to fully chaotic. Their quantum counterparts attracted 
much interest starting in the 80's, both from the point of view of random matrix theory and the 
semiclassical periodic orbit theory. In the latter approach one uses trace formulae of the type first 
derived by Gutzwiller Q and Balian and Bloch |2|, |3[. 

During the last two years, following the route opened by Ref. Q, a number of studies (see Refs. 
[^, ^, 1^, IH) have concentrated on additional contributions to the trace formula linked to diffractive 
effects near regions where the classical Hamiltonian flow is discontinuous. These zones of discontinuity 
are known as "optical boundaries" in the literature. They lead to contributions from non-classical 
(so-called diffractive) orbits hitting a corner of the billiard or creeping around a smooth boundary. 

Apart from the noticeable exception of Ref. all the work quoted above is based on Keller's 
"geometrical theory of diffraction" (GTD, see e.g. |^), i.e. on an extension of geometrical optics 
which accounts for diffractive effects. Keller's approach fails when the diffractive trajectory is very 
close to an optical boundary, or equivalently when the diffractive orbit is close to become an allowed 
classical trajectory (this will be clarified in the text of the paper). In the present work we use a 
uniform approximation for the Green function which does not have this drawback. This allows to 
derive relatively simple formulae which are uniformly valid. The method is applied to billiards whose 
boundary has a slope discontinuity and thus we restrict our study to wedge diffraction effects. To our 
knowledge there does not yet exist a uniformly valid formula for the contributions of creeping orbits 
(despite the progress made in Ref. 0). 

The theory of uniform approximations for wedge diffraction has a long history which begins with 
a famous paper by Pauli |10|. In the late 60's and in the 70's the problem has been studied in detail. 
Much literature has been devoted to several types of approaches remedying the deficiency of the 
geometrical theory of diffraction. The approach most widely used is known as "uniform asymptotic 
theory" and was developed by Ahluwalia, Boersma, Lewis and coworkers in Refs. |11, |l^, |l3|, ^]. 
We have chosen here a technique more closely related to the original work of Sommerfeld and Pauli. 
It relies on an extension of the method of steepest descent due to Pauli which was carefully studied 
on a general setting by Clemmow |15|. The method, due to Kouyoumjian and Pathak, is known 
as "uniform theory of diffraction" and is exposed in Refs. and [17|. Note that we apply the 
uniform approximation only to orbits with a single diffractive point. The treatment of multiple wedge 
diffraction is increasingly more involved as can be seen in work on double diffraction by half-planes 



(see Refs. [18, 19, ^) or wedges |21|. To our knowledge there does not exist to date a general uniform 
approximation for multiple wedge diffraction. 

The paper is organized as follows. In the next section we recall the exact solution of the infinite 
wedge problem, derive a uniform approximation for the Green function and compare it with the result 
obtained from GTD. In Sec. 3 we use the Green function obtained previously to derive contributions 
to the trace formula which are uniformly valid. Readers mostly interested in the final result can skip 
this part and go directly to Sec. 4, where we discuss the previously obtained formula and several of its 
limits. In particular, we show that this formula has the appealing feature of interpolating between the 
semiclassical results of periodic orbit theory and the formulae obtained in Refs. ^, ^. Sec. 5 contains 
numerical applications for several simple billiard systems. In some cases GTD gives reasonable results, 
but in other cases the uniform approximation has to be used in order to describe the Fourier transform 
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of the spectral density correctly. Finally we discuss our results and possible extensions in Sec. 6. 



2 The Green function of an infinite wedge 

In this section we consider an infinite wedge of interior angle 7 (7 G]0,27r]) with Dirichlet boundary 
conditions and derive several approximations for the Green function. 

2.1 The exact result 

The exact solution of the problem was first given by Sommerfeld for a wedge with 7 = 27r (a half line) 
and an incident plane wave, see p^. The solution of the general problem is easily inferred from his 



approach, a complete treatment is given for instance by Carslaw in Refs. |23, 24 1. We recall here for 
completeness some properties of the solution. 

The Green function G.y(r, r ', E) of the problem in dimensionless units is a solution of: 



iAr- + E)G-,if,f',E) = 5{r-r'), (1) 
G-y = if f or r ' are on the boundary . 

Choosing a system of coordinates with the origin at the vertex and the polar axis along one of 
the boundaries such that 9 and 9' are in [0,7] (see Fig. (2a)), one can write the following integral 
representation for the exact solution: 



G^(r, r", E) = g^{r, r', 9' - 9) - g^{r, r', 6' + 6) , (2) 



with 



8ttN Ja+b 



H^^^ {k\lr^ + r''^ - 2rr'cosz^ 



1 — exp 



z 
-i- 



N 



(3) 



In (P) and in the following the angles 9 and 9' always appear in the combination 9' zk9 and we will 
denote = 9' — a 9 (fx = ±1). Other quantities appearing in (|3[) are = 7/vr, k = ^/E which is the 
modulus of the wave- vector, Hq^ which is the Hankel function of the first kind (see ||25|), and A and 



B which are the contours in the complex plane drawn in Fig. 1. In this figure one can further see the 
poles of the integrand corresponding to {z — (pa-) /N = 2mr (with n £ TL) - they appear as black points 
- and branch cuts linked to the square root argument of The shaded areas are zones where 

the integrand increases without limit when one goes away from the real axis (this is easily checked 
by using the leading asymptotic term (^) of the Hankel function). The integration contour is quite 
arbitrary as long as it goes to infinity in the indicated unshaded regions. 

Essentially, is a superposition of free Green functions (with complex angles z). Considered as 
a function of (p^, it has periodicity 2-11 N = 27 and this ensures that the Green function @) satisfies 
the boundary conditions for 9 = and 7. By moving the contours A and B towards the real line 
z G [— TT, vr] and taking into account the poles of the integrand one obtains 
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Figure 1: Integration contour in the complex plane for formula (P). The shaded areas are zones where 
the integrand diverges when going away from the real axis. The black points are poles and branch 
points of the integrand. The thick lines are branch cuts. 



(7^(r,r',0^) = J2' ^0^^ ( fc^r^ + r'^ - 2rr' cos{<j)^ - 2n-f) ) + /i^(r, r', 0^) , (4) 



4 „ 



where after a change of variable can be written in the form 



cos [ — — — j — cos(7r/A') 

The first term on the r.h.s. of (Q) contains the contributions of those poles of the integrand of (^) 
which lie between — vr and tt; the prime indicates that the summation is restricted to values of n such 
that — vr < (j)cr — 2n7 < vr. If (p^ is exactly equal to itvr + 2nj then the corresponding contribution to 
the summation has to be divided by 2. In ^ the contour can be modified as long as no pole of the 
integrand is crossed. A further requirement is that the part of the contour extending to infinity has 
to start at —ioo with a real part in [0, vr[ and to extend to ioo with a real part in ] — vr, 0]. 

The discrete summation in can be interpreted as arising from allowed classical trajectories. For 
instance the case (pj^ = 9' — 9 and n = gives a contribution —{i/ 'i)H^\k\r — r'|) and corresponds 
to the free propagation from r ' to r. The other terms in the summation correspond to trajectories 
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experiencing specular reflections on the boundaries (this is illustrated in Fig. (2b)). If a = 1 (resp. 
a = —1) the orbit has an even (resp. odd) number of reflections. These orbits correspond to successive 
applications of the method of images and their contribution is known as the geometrical term in the 
literature. When inserted back in Eq. (Q) they give a term which will be denoted Ggeo{r, r E) in the 
following. 

If the angle 7 is of the form ir/p [p G IN*) then sin(7r/A^) = and the term /i^(r, r', (/iq-) is zero: 
the geometrical term alone is enough to fulflll the boundary conditions. This is due to the fact that 
in this case the Green function can be determined by the method of images. If 7 7^ ^r/p, then 
corresponds to the contribution from diffraction. Hence the total Green function can be written as a 
sum of a geometrical plus a diffractive term: 



G^{r,r',E) = Gg^r^r' ,E) + Gdiff{?,?' ,E) , (6) 
with Gdiffif,f',E) = h^{r,r',e' -e)-h^{r,r',6' + e) . 

2.2 Geometrical theory of diffraction 

We derive now a simple approximation for Gdiff- We first replace the Hankel function by its asymptotic 
form for large arguments (see [p5|): 



i7{i)(2) ~ W4 ^hen |z| » 1 . (7) 

V TTZ 

The same approximation is used also in all the following for the geometrical and diffractive Green 
functions, i.e. for all the terms of Eqs. (^) and (^), the assumption being that the distances measured 
along the paths (classical or diffractive) going from r ' to r are large compared to the wavelength 
A = 211 /k. Then in the integral defining /i^ there is a saddle point of the exponent at z = and a 
steepest descent approximation yields: 

_i_ sinin/N) gifc(r + r') + i7r/2 

h^{r, r , 0,) « cos(0,/iV) - cos(7r/iV) ' 

Incorporating this result into the expression (|6|) for G^iff one obtains a formula which can be cast 
into the form: 

Gdiff {r, f', E) ^ Gscif, fo, E)V{9, e')Gsc{fo,r' , E) , (9) 

with 



V{9,e') 



2 TT 



vr 



+ 



-1 



-1 



cos ■ 



N 



IT 



COS ■ 



(10) 



4 sin(7r/iV)sin((9/iV)sin(6'7A^) 

(cos — — COS — — )(C0S — — COS — — - 

^ AT N N N 
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In dll) Gsc is the free Green function evaluated using and ro is the point at the vertex (see Fig. 
(2a)). Expressions (|9|) and (p!o| ) give the diffractive part of the Green function in the "geometrical 
theory of diffraction" (see [y). They have the simple interpretation as being the contribution of a 
(non-classical) diffractive trajectory going from r ' to tq and then from rg to r (see Fig. (2b)). Using 
this approximation one can derive a trace formula for the spectral density which accounts for diffractive 
effects in the GTD approximation (see Refs. ^). 

( a ) ( b ) 




Figure 2: Figure (a) displays the notations used in the text. Figure (b) shows two classical trajectories 
(solid lines) and the diffractive orbit going from r ' to r (dashed line). 

The quantity T>{6, 9') is known as the diffraction coefficient. It is zero if 9 (or 9') is equal to or 7, 
or if 7r/7 is an integer. It diverges on an optical boundary, i. e. if r and r ' are such that the diffractive 
orbit is the limit of a classical trajectory. This is illustrated by the simple case of diffraction by a sharp 
wedge (7 > 37r/2) in Fig. 3. In this case there are two optical boundaries represented by dashed lines. 
They correspond to 6* = 6*' + vr (i. e. = 6*' - 6* = -vr) and 9 = tt - 9' {i.e. (p- = 9' + 9 = tt). In the 
terminology of geometrical optics, the first optical boundary separates the illuminated and shadowed 
regions for direct rays when the point r ' is considered as a light source (boundary between regions 
II and III), and the second optical boundary separates the illuminated and shadowed regions for rays 
that are reflected on one side of the wedge (boundary between regions I and II). If flies near one 
of the optical boundaries then the diffractive path is almost an allowed classical trajectory, and if f 
is moved onto an optical boundary then the diffractive path coincides in this limit with an allowed 
classical trajectory. 

Looking in more detail at the origin of the divergence, one sees from Eq. ^ that it occurs when 
there exits an integer n such that = ±vr + 2nNTT. In this case there is a pole z = (pu^Ti — 2nNTT = 
in the integral representation of the diffractive part (^) which is at the same position as the saddle 
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region III 



Figure 3: Optical boundaries (dashed lines) for an initial point f in the case of diffraction by a sharp 
wedge (7 > 37r/2). The transition regions (for r' = lOA) around the optical boundaries are shaded. 



point 2 = and thus the saddle point approximation breaks down. More generally, the geometrical 
theory of diffraction is only valid if all poles are sufficiently far away from the saddle point z = 0. This 
can be interpreted in terms of the physical trajectories in the system because, in Sommerfeld's solution 
^^j, the saddle point corresponds to the diffractive orbit and the poles correspond to geometrical 
orbits. 



2.3 A uniform approximation 

As seen above, one has to refine the steepest descent evaluation of in case that there is a pole of 



the integrand near the saddle point z = 0. This was first done by Pauli [10| and we present here a 
slight modification of the original procedure |T^. In a first step one can separate poles which are 
possibly near one another by using the identity: 

2sin(7r/iV) _ 1 1 ^^^^ 



<^os{^^ -cosin/N) tan( — j tan ( — 
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Hence hy in Eq. (||) can be rewritten as 

h-y{r', r, 4>a) = u^,+ {r', r, (j)^) - u^-{r', r, (j)^) , (12) 



where 



X r+too H^^^ ky + r'^ + 2rr' cos z 



and 7/ = ±1 is a new index. 

If one denotes by ng-,,; the nearest integer to (c/^o- + rj 7r)/{2^) then z = —{(pa + r/ vr) + 2 no-,r)7 is the 



pole of the integrand of ( [l^ ) which is nearest to the saddle point z = 0. Thanks to the separation (12) 
the next pole in the integrand of ([l3[) is at distance 2j and its effect can safely be neglected if 7 is not 
a small angle (this will be assumed in the following). According to the method of Pauli one rewrites 
the integrand by multiplying numerator and denominator by a function imitating the behaviour of the 
original denominator but in which the z and (pa parts are separated. This procedure is not unique, 
it corresponds to a specific choice of a uniform approximation as will be discussed below. The choice 
for the function is rj -v/2sin(z/2) + Ua^r], where Ua^r) is a measure of the separation between the saddle 
point z = and the nearest optical boundary: 

aa,ri = \/2 cos(% - ria^nl) with n,^ „ = nint [^^-^^-^] G ^ • (14) 
2 27 

Using the asymptotic formula (0) for the Hankel function one obtains: 



p-i7r/4 „+ioo r,ik\l r"^ + r'"^ + 2rr' (zos z 

u,,,{r, r', (Pa) - -j= / . dz ^ F^'^z) , (15) 

8jv2TTk J -too rjv 2sm.[z/ 2) + aa^r^ 

where F"'^{z) is a smooth function at z = 0, even in the vicinity of an optical boundary (when 
aa,r) 0, see (|2^)): 

pa,r, ^^-^^ r/ ^/2 sin(z/2) + a^,^ 



2 + /2 ^ 2rr' cos zf'^ tan {^-±^^±^^ 



[r 



Note that the integrands of (|T^) and (13) both have the pole z = —{(pa + Tjir) + 2na^nl next to 
the origin as mentioned above. 

Now ( p^ ) is evaluated along the steepest descent path at z = by a change of variable z = 
exp(i37r/4) with t G IR (the factor riy/2 is here for convenience). The smooth, non-singular part 
F'^''^ of the integrand is simply evaluated at i = 0, and the phase of the exponential function and the 
denominator are expanded in the vicinity of the origin: 

r]V2sm{z/2) + aa,r^ = aa,^ + t e^^^/^ + 0{t^) , (16) 

and 
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ik\l + r'"^ + 2rr' cos z = ik{r + r') - , f + 0{t^) . (17) 
Hence approximated by: 

Q,ik{r + r') - m / A „oo exp -krr't^ / {r + r') 

u^,r,{r,r',<P^) ^ . / dt — . (18) 

8j^7Tk{r + r') tan(^^^t^) t - a.,,eW4 



After the expansion (jl^) of the denominator, the pole in ( p^ is only approximately equal to 
the nearest pole in (^), but they coincide when the pole approaches the stationary point t = 0. 
As noted above the choice of the uniform approximation which leads to Eq. (|l^) is not unique (as 
discussed by Clemmow, who calls it a "partial asymptotic expansion", see [^). For example, another 
choice of a uniform approximation can be obtained by making a change of variable which transforms 
the exponent in (|l^) such that it becomes an exact quadratic function, and multiply denominator 
and integrand by a function which is linear in the new variable. Then one obtains Eq. (|l^ ) with 
a different definition of Uf^^ri which is expressed in terms of a "detour parameter" such as used in 
the "uniform approximation theory" (see Refs. 12, 13, Q). In practical applications, however. 



the differences between different uniform approximations are small. We would like to add that the 
uniform approximations can be further improved by including also sub-leading terms of the asymptotic 
expansion of the Hankel function and also higher order terms of the expansion of the integrand (see 
Refs. 11, H]). However, numerical checks show (see Fig. 4 and below) that such refinements are not 



necessary here. 

We continue now with the integral ([l^) which can be recognized as an integral representation of the 
modified Fresnel function K (see Appendix A and Eq. (|A5| )) and the final expression for the uniform 
approximation for u^^^ is: 

^ Qik{r + r')+i7r/A , , / rk^\ 

This expression remains finite on the optical boundary (p^j = —r] tt+2 rifj^r)!- As an optical boundary 
is crossed, ag-^r; goes through zero and changes sign and one has: 

I ^(T,ri I 



, , Pa+rj-K 

tan 



i r]V2N sign (a^,?;) when aa,r] . (20) 

2N ~) 

Hence although the problem of divergence has been eliminated one arrives at a final form which is 
discontinuous. This was expected: the exact terms ^ and ( |l3[) already have this behaviour; because 
of the separation (^) of the total Green function into a geometrical and a diffractive term, each 
contribution (Ggeo and Gdiff) is discontinuous at the optical boundary, but their sum is continuous. 

As a resume of the results of this section we write down the uniform approximation for the diffrac- 
tive part of the Green function which is a sum of four contributions: 
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1 pifc(r + r') + m /A 

G,.ff{r,r\E) . E "A'^^I^x ^ I 1 > (21) 

'7r/i;(r + r') a,r,=±i tan 

where 4)^ = 0' — cr {6 and 6*' being chosen in [0,7]) and aa,r) is defined in ([l^ . 

In the remaining part of this section we present some numerical results illustrating the accuracy 
of the uniform approximation and a failure of the GTD approximation. If the next optical boundary 
is sufficiently far away one can replace the modified Presnel function in (21) by the first term of 



its asymptotic expansion (A4) and this leads to the GTD result (p|,[lO|). Roughly speaking, this 



approximation is good when the argument of the ET-function is greater than 3 and it fails when 
the argument is less than 1.5. This puts a limit on the use of the geometrical theory of diffraction 
illustrated in Fig. 3: inside the dashed areas around the optical boundaries one has to use the uniform 
approximation (these zones are known as "transition region" in the literature). The figure has been 
drawn for the case r' = 10 A (A = 2ir/k), and the transition regions are larger if one goes to smaller 
values of r'/A. In the limit r ^ r' ^ A, the transition width around an optical boundary at distance 
r from the apex is proportional to r^fXpP (relying on the weaker assumption that r, r' ^ A one can 
show that it is proportional to [A(r^ + rr')/r']^/^). Outside of the transition region expression (p|JTO|) 
is valid and shows that G^i// is a small correction to Ggeo- But near the optical boundary the two 
terms are of the same order and exactly on the boundary the two discontinuous contributions to Ggeo 
and Gdiff have exactly the same amplitude. 

The comparison between the uniform approximation (^l|), the geometrical theory of diffraction 
(p|JTO|) and the exact result (^,^) for Gdiff is made quantitative in Fig. 4. In this figure one considers a 
wedge of interior angle 7 = 110°. The source point r ' is fixed at 0' = 60° and r' = 5A. The observation 
point r is at fixed distance from the vertex (r = r') and Q scans the interval [0, 7]. The modulus of G^if^ 
is then plotted as a function of ^. In the figure one cannot distinguish the uniform approximation from 
the exact result. The geometrical theory of diffraction diverges on the optical boundaries (represented 
as dashed lines in the upper part of Fig. 4). Furthermore it is in clear disagreement with the exact 
result for all values of Q. Hence one can infer that a trace formula based on Eq. (^) will not correctly 
describe the spectrum in cases such as presented in Figure 4. 



3 Diffractive orbits in the trace formula 

We consider now a closed two-dimensional region B with a boundary dB smooth everywhere except 
at a finite number of points where its slope is discontinuous. The spectral density d(A;) of this system 
has semiclassical contributions from periodic orbits as well as from diffractive orbits. The latter ones 
are closed orbits which have a finite number of points on vertices of the billiard (we call these points 
diffractive or corner points in the following) and follow the law of geometrical optics between two 
diffractive points. Within the framework of the geometrical theory of diffraction, the contribution of 
a diffractive orbit ^ to the level density has been derived in 1^, ^ and it is given by 



n 



\ y/8TTk{Mi) 



12 



cosikL u p) 

\ 2 4 



(22) 
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Figure 4: Modulus of Gdiff{r,r' ,E) for fixed r' and r in a wedge with 7 = 110° (r = r' = 5A and 
9' = 60°). ^ scans tlie interval [0,7]. The upper part of the figure displays the geometry considered, 
the optical boundaries appearing as dashed lines. In the lower part, the solid line is the exact result 
(^Q, the long dashed line is the uniform approximation ( ^T|) and the short dashed line is the GTD 
result (||,|0|). 
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Here Lq and L are the primitive and total length of the trajectory, respectively, p is the number 
of diffractive points, Vi is the diffraction coefficient in the i-th. corner (cf. Eq. (pX|)), (Mj)i2 is the 
(12)-element of the stability matrix at unit energy for a part of the trajectory between two corners, 
and v is the number of conjugate points plus twice the number of reflections on the boundary between 
corners. 

According to ( p2[) each corner point decreases the contribution of a diffractive orbit by an order 
©(fc^^/^). This is correct only if the diffractive trajectory is sufficiently far away from the optical 
boundaries in every corner point. In the opposite case that the trajectory lies on an optical boundary 
in every corner point it can be shown that its contribution is of the same order in k as that of a regular 
periodic orbit (see below). In the following we will go beyond the GTD approximation and derive a 
contribution to d{k) from diffractive orbits with one point in a corner that interpolates between these 
two regimes. We will use for this purpose a method of uniform approximation similar to that exposed 
in the previous section. 

The starting point of our derivation is the boundary element method. It is a reformulation of the 
quantum mechanical eigenvalue problem in terms of a Fredholm equation of the second kind for the 



normal derivative of the wave function on the boundary (see e. g. Refs. ||26|, ^ |28|, ^] for discussion 
and application in the context of the trace formula). More specifically if we denote by f{s) a point of 
the boundary with curvilinear abscissa s and by u{s) the normal derivative of the wave function at 
this point, one has the following integral equation for the case of Dirichlet boundary conditions: 

u{s') = -2[ dsu{s) dn'Go{r,r',E) , (23) 
JdB 

where r = r{s), r' = f'{s'), h' is the outward normal vector to dB at point f', d^i is the projection of 
the gradient onto n', and Go{f,f',E) = —{i/4:)HQ'\k\f—f'\) is the free Green function. The integral 
relation ( p3|) has non-vanishing solutions u{s) only if 



det (/-Q(A;)) =0, (24) 
where / is the identity and Q{k) is an integral operator which, when applied to the function u{s), 



gives the r.h.s. of Eq. (23). The zeros of (|24| ) are the exact quantum energies of the system and the 



oscillatory part d{k) of the level density can be expressed as 

-17 1 ^ 1 W 

J(A;) = — lndet(/-Q(fc)) = -9 V- — fTrQ"(fc)l , (25) 

n=l 

with 

TTQ^{k) = {-2r I dsi..Jsndn,Go{r2,ri,E)dn,Go{r3,f2,E)...dn„Go{n,rn,E) . (26) 
JdB 

For a system with a boundary B that is smooth everywhere, the integrals in ( p6[ ) can be evaluated 
in stationary phase approximation. In this way Tr Q"^ is expressed in terms of a sum of contributions 
arising from periodic orbits with n specular reflections. Inserting this approximation into Eq. ([25|) 



yields Gutzwiller's trace formula as shown, for example, in Appendix B. 

The standard approach described above is not convenient for deriving contributions of diffractive 
orbits, since the diffractive effects of the corners are hidden in this formulation. Instead, we use a 
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modification of the boundary element method by formulating it in terms of a Green function accounting 
for the diffractive effects of a corner. 

We restrict to the consideration of diffractive orbits with a single corner point which are not 
influenced by the other corners of the billiard. To obtain the contribution of such an orbit to the trace 
formula it suffices to include the diffractive effect of only one corner. We further restrict in this section 
to corners in which the limit of the curvature of the boundary is zero when the corner is approached 
from either side. Modifications caused by non-vanishing curvature are discussed in Appendix D. 

Let us first consider a simple billiard system which is bounded by a wedge of angle 7 and an 
additional smooth curve C which connects the two sides of the wedge (such as represented in Fig. 6 
for instance). One can derive an integral equation in terms of the Green function of the infinite 
wedge - analogous to equation (23) - in which the integration is restricted to the curve C. The 
oscillatory part of the spectral density is then again given by (25) where the trace of Q"' now has the 
form 



TrQ"(A;) = (-2)" / dsi . . . dsnd^, G^{f2,fi, E) d^, G^{r^,r2^E) . . . d^^ G^{n,rn,E) . 
J c 



(27) 



Note that similar techniques have been used in Ref. Q for deriving diffractive contributions in the 



Sinai billiard and in Ref. |31] for reformulating Fredholm's theory in the case of triangles. 

The Green function can be split into a geometrical and a diffractive part as has been done in 
section 2 



G^(f, f', E) = Ggeo{r, r', E) + Gd^ff{r, f, E) , (28) 

where the diffractive Green function is given by expressions (|,|1|,|13D. Inserting Eq. (||) into d^) 
results in 2"" integrals. The stationary points of these integrals correspond to periodic and diffractive 
orbits of the billiard system (and possibly also to ghost orbits as in the case of billiards without 
corners), and the number of points in a corner of a diffractive orbit is determined by the number of 
diffractive parts G^iff appearing in the integral. Since we restrict to orbits with one point in a corner 
we can replace {n — 1) of the Green functions in (27) by their geometrical part. This can be done in 
n ways which cancels the factor 1/n in (|25|). Then the contribution to the level density from orbits 
with n refiections on the boundary C and one point in a corner are contained in 



^ (—2)" d f 

dY'{k) = —^-TT I dsi... dsn «9fti G^{r2,fi,E)dn^ Ggeo{n, ^2, E) . . . df,^ Ggeo{ri,rn, E) , (29) 

vr ak jc 

an example is given in Fig. 5. 

Let us discuss Eq. (^) in more detail. The contribution of a diffractive orbit is obtained by 
evaluating the integrals in the vicinity of the stationary points, i. e. in the vicinity of the points of 
specular reflection of the orbit on the boundary. If one approximates G-y in the framework of GTD this 
results in Eq. (p2|) (with p = 1) for the contribution of the diffractive orbit. In the following we will 
improve on this method by using a uniform approximation for the Green function G-y. In both cases 
however, only local information about the reflection points and the corner enters the approximation. 
It is then obvious how the expression ( p9| ) has to be modifled in order to derive the semiclassical (or 
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Figure 5: Typical path contributing to (p9|). A precise definition of the angles Oj can be found in 
Appendix B. 

uniform) contributions to d{k) for more complicated diffractive orbits in billiards with several corners: 
for every straight part between two reflection points a free Green function has to be included, and 
for every part of the trajectory which hits a corner between two reflections a Green function for an 
infinite wedge with the same angle. The reason why Ggeo is appearing in Eq. (^) and not Go is that 
in the above formulation we consider only reflections on the part C of the boundary, and Ggeo takes 
care of the reflections on the wedge part. Hence the total number of specular reflections on dB in Eq. 
( p9[ ) may be greater than n. 

We continue now with the further evaluation of Eq. (|2^) which we perform in the case n 7^ 1. The 
calculations for n = 1 can be treated by identical methods and yield the same final result. 

In ([29|) , (n — 2) boundary integrals can be evaluated by applying the composition law ( ^l]) for 
Green functions which is derived in Appendix B: 

dA,Gscifl,r2, E) ^ (-2)("-2) I ds^... dSn dn, G geoifs , r2 , E) ...dn^ Ggeoifx^Tn, E) , (30) 

and consequently 
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df\k) = ^^^^ I dsi dS2 dn, G^{f2,ri,E)dn, Gsc{fl,r2,E) . (31) 
vr ak Jc 

Here Gsc is the contribution to the semiclassical Green function from trajectories with (n — 2) 
reflections on the boundary curve C (and possibly further reflections on the wedge part of the boun- 
dary) : 

Gsc{ri,r2,E) = ^ ^==== exp{ifc/ - i-P - i— } , (32) 
^ V87rA;|mi2| ^ 4 



where m is the stability matrix (see Appendix B| ) and / the length of the classical orbit going from 



r2 to r*!. i) is the number of conjugate points plus twice the number of specular reflections on the 
boundary. We use here u and lower case letters for ra and / in order to distinguish these quantities 
from those of the whole diffractive orbit. The normal derivative of the Green function is given in 
leading order by 

dn2Gsc{ri,r2,E) 7^ ik cos a2Gsc{r I, f 2, E) , (33) 

where 02 is the outgoing reflection angle at r2 (see Fig. 5). 

In the following, we will consider the contributions of the geometrical and diffractive parts Ggeo 
and Gdiff to the Green function G^ in Eq. (^) separately. As discussed above, the geometrical part 
will yield the contributions of periodic orbits. The reason why it has to be included also for the 
derivation of the contributions of diffractive orbits is that both Gdi / / and Ggeo are discontinuous at 
the optical boundary (see the discussion in Sec. 2). For that reason the boundary contribution of Ggeo 
which arises from this discontinuity has to be included in order to cancel the analogous contribution 
of Gdiff. 

3.1 The diffractive contribution 

Prom (^,|7, 1^,13) the diffractive part of the Green function G^ can be approximated by 



/ 2 e~^^/^ f'-°° expji/cwr^ + r| — 2rir2 cos 2} 

Gdiff{f2,fi,E)^ J2 ^m^—Tr / / , , , ^ A^^) 

.,.=±1 Vvrfc I67 7-.^ (,2 ^ ^2 _ 2r,r2 cos z)V4 tan + + 



2N 

where ri (resp. r2) is the distance from fi (resp. r2) to the diffractive point, and (p^j = 61 — a92 (see 
Fig. 5). Similarly to (|33|), the normal derivative yields a factor ik cos ai. We insert ( |3^ and one 
contribution from ( |32[ ) into (^Tj) and consider the contribution from the vicinity of a stationary point 
which is chosen as origin of the s- variables. The main contribution to the ^-integral comes from values 
near z = and the exponent is expanded in z up to second order 



di4iff{k) 



-2f 



vr 



d X - 

■Tk ^^/^^ 



;37r 



a,ri= 



nk I67 



\^oTTk\rn^\ 



{ik)"^ cos ai cos a2 
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«oo exp{iA;(/(si, S2) + r-i(si) + r2(s2) - 2?^^^^^^^)} 

dsi ds2 / dz 7 — —7 — ^ — , (35) 

Vri + rs tan ( — 

where the index ^ labels the difFractive orbit. A stationary phase approximation of all integrals would 
yield the contribution of the diffractive orbit in the GTD approximation. This approximation diverges 
at an optical boundary. In order to obtain a finite uniform approximation the effect of the nearest 
pole to 2; = has to be included. We treat this pole again by the method of Pauli 



.^^+^^sin^±^ 
1±A^+JE\ ~ tan { l+^+iZZl') a,,„ + r?^/2 sin 



tan ^ ^ ' j tan ^ " ' j a^,^ + r/V^ sm — ^ 

1 a„ 



tan^-^j^ ) ''-'V + V ^ 

where (/>o-,o is the value of 4>a at the stationary point, Acft^j = cj)^ — (p^ Q and aa^n is evaluated using ([1 
at the stationary point. Inserting p^) into (53) we obtain 



(36) 



d ^-^ fc cos ai cos q;2 exp{— i^i/} 



dk^^^^, 87r27V(ri + r2)|mi2| ( 



too exp{iA;(/(si,S2) +n(si) +7-2(52) - tTZ^-^tVt^^)} 
dsids2 I dz /^"^^ + '^^> . (37) 



—00 J —lOO 



aa,r, + viz + A</>a)/\/2 

The quantities l{si,S2), ri{si) and r2(s2) are now expanded up to second order in si and S2- The 



expansion coefficients can be obtained from (B8). Furthermore we expand A(p(j up to first order in si 



and S2- A<pcr{si, S2) ~ sicosai/ri — cJS2 cos a2/r2. After a substitution 

V2 V2 



si -r] si , S2 -r] S2 , z -rjV^z , (38) 

cos ai cos a2 



we obtain the following expression 



d /cexp{— z^P} acr,n 



dk ^ 7r27V8(ri+r2)|mi2| , ( <Pofi + W 

2N 

00 MOO exp{iA;(/ + ri + r2 + a' si + h' + j^siS2 - cz^)} 
dsi ds2 I dz , si S2~~ 

(39) 



00 J —lOO 



7 -I- i-i — — n 

^ ' Tl T2 (^1^ 



where 

a = ^ + — , h = + — , c = . (40) 

mi2 /ticosai ri mi2 /i2Cosa2 ^2 n + r2 
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Here the quantities ri, r2 and / without argument denote the values at the stationary point. The 
derivative with respect to k in ( p9| ) yields in leading order a factor iL where L = ri +r2 + / is the length 
of the diffractive orbit. In the next step we simplify the integrals by applying a transformation of the 
s variables such that the denominator of the integrand depends only on one of the new s variables 

s = — — a— , s' = dr2Si + (1 — a d) riS2 ■ (41) 

The form of s' is chosen such that the Jacobian of the transformation is one and the value of d is 
determined by the requirement that the exponent in the integrand has no mixed quadratic term s's. 
The evaluations are done with Maple and result in 

d!i4iff{k) ^ -?fi ar] 



2 I UJ2 „^2\ 



f°° , , , r~ , exp{ik(as^ + bs'^ - cz^)} 

/ ^^^^ / . TTT^n ' 

J —oo J —loo ^ -r o "(T.rj 



where 



Mi2C 



Mi2 - c(Tr M - 0-2) 



d 



Mi2 — c(Tr M — a2) ' mi2rir2C ' 

c{ariRi cos aimu — 2arimi2 + crRi cos aimi2 + r2Ri cos ai) 



(43) 



r2i?i cos ai [Mi2 - c(Tr M - a2)] 

M is the stability matrix of the diffractive orbit, i. e. the stability matrix of the classical trajectory 
starting from the corner going through r2 and fi and back to the corner. There is a relation between 
D, the signs of a and b and the Maslov index i/ of the diffractive orbit 



expj-i-z^ + ^-c^a + ^^o-fej = i ex.p{-i-u} , 



(44) 



where aa = sign(a), ab = sign(6). Hence v is equal to P plus the number of negative signs of a and 



b (modulo 4). This can be seen, for example, by evaluating the si- and S2-integrals in Eq. (3£) by 
the stationary phase method. Then from the composition law of Green functions the Maslov index 
u of the whole diffractive orbit is obtained by successive applications of ( B12| ). Since a stationary 
phase approximation after the transformation ( pT[ ) has to yield the same result, relation ( |4^ ) follows 
immediately. 

The integral over s' can now be evaluated, and the double integral over s and z is calculated in 
[Appendix C| (Eq. (^). The result is 



Lexp{ikL — i^v} 



I 



a,^±i " lV8{ri+r2)\mi2b{a-c)\ ( 0^.0 + 



ii7r(l+(Ta+cr6+TCT)/4 



ikacal 
expj 



a — c 



erfc{ I aao-,r; I 



i[a — c 



^ } -erfc{|a^,^| 



^ a — c 



(45) 



17 



where To- = sign(ac/(a — c)) = sigii(Mi2/(Tr M — o"2)). As will be seen in the following section, the 
first error function in ( p5| ) is the contribution from the discontinuity of Gdiff which is cancelled by 
the corresponding contribution from Ggeo- 

3.2 The geometrical contribution 

For a given a and r] the geometrical orbit that corresponds to the nearest pole arises when (j)u—2na,'q^ = 
—r]TT and it exists if r]{2na,r]1 ~ 4'a) < tt- This can be re-expressed in the form 

aa,n = v2cos( —) > -V2sm — - — « . (46) 

2 2 y 2 

In the following we derive the contribution from the discontinuity of Ggeo- For that purpose we 
apply exactly the same approximations to Ggeo that were used for Gdi / / • This is done by writing Ggeo 
in the form 



Ggeo 



a,r]=±l 
o-,77=±l 



2 e^^/ 4 exp{iA; J r\ + T\- 2rir2 cos((?!>^ - 2n„^^-f)] 



(rf + rl 



2r\r2 cos( 



(47) 



exp{z/i; w rf + r2 — 2rir2 cos z} 



irk I67 / 



(r 



2 + ^2 _ 2r^r2 cos z) V4 tan (^£±i^lM^ ' 



where A = aa,n + \/2 sin(?7A(/>o-/2). The integration contour of the z-integral encircles the nearest 
pole to z = counter-clockwise. The expression (47) differs from Eq. ( |3^ ) only by a factor (—?]), 
the 0-function and the integration contour. We repeat now all the steps from Eq. ( |35| ) to (p^. The 
only difference is an multiplicative factor {—rj) which results from the substitution z —r\\f2z (it 
didn't appear previously because of the different integration contour). We arrive at an expression 



corresponding to Eq. (42) 



ar] 



kLcxp{ikL — i^D} aa,Ti 



° ds ds' idz e(a,, - , fM^kias^ + bs--cz^)} 

yo J Z -\- S CLcr^TI 

The triple integral is denoted by /. The integrals over s' and z can now be evaluated and result in 



2Tri 



ds expjifeas^ — ikc{s — a^^r])'^} 

I TTZ etc 

2^m/ 7T / ds exp{ik(a — c)(s -\ 

V kb J -00 a 



kb J —00 
rO 



(49) 
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We are interested only in the boundary contribution of the geometrical part. The expression (|4^ ) 
contains in general also contributions from periodic orbits. This is the case if the integration range in 
( ^9|) contains a stationary point, i. e. if aa^^r]/ {o- — c) is positive. Then the stationary point contribution 
has to be subtracted which corresponds to a subtraction of the integral from — oo to oo. We obtain 
for the boundary contribution 



I' = — sign(ao-j7) To-27riW — / ds exp{ifc(a — c)(s + 

V f^t) Jo 



aa, 



a — c 



a — c 



-sign(a^_^) Tc 



k V|6(a-c)| 



1 Qi^i'^+"'^+^b+^'^) 



a — c 



}erfc{|aao-,, 



i{a — c) 



} , (50) 



where Tq- = sign(ac/(a — c)), as before. Substituting /' for the triple integral in (^) yields 



lV8{n + r2)\mi2b{a - c)\ ( (pa,o + W 

2N 



a,ri=±i /v"Vi ' 'zjr'-'Lz^y^ witan 



a — c V t{a — c) 



Comparison with Eq. ( ^5| ) shows that this contribution exactly cancels the first error function in 



3.3 The joint contribution 

We now can write down the final formula of this section. By using the definitions of a, b and c the 
sum of ( |45| ) and (^) can be written in the form 



L expjifcL — i^fj-a} 



''cr,ri I 



-5R TIT 

..=±1 '^^ V|TrM-a2| 2Ar^tanf^^ 



exp{- 



ikalj^Mu 
Tr M - cj2 



} erfc{|a, 



kMi2 



■a,ri I 



i(TrM-(T2) 



} 



(52) 



where we dropped the second index of 0o-,o for simplicity of notation. Furthermore, /io- = v-\-{1—(t)+Ko 

M12 



= 1 — 2Ka, and Kq- is defined as 



if 



1 if 



Tr M - cj2 



M12 



> 



< 



(53) 



Tr M - cj2 

We recall that M in (52,^) is the stability matrix (at unit energy) of the classical trajectory 
starting and ending at the corner point, v is the number of conjugate points of this trajectory, plus 
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2 times the number of specular reflections. The definition of in terms of u and Ko- is similar to 
the definition of the Maslov index of a periodic orbit in terms of that of the Green function |33|. The 



cj-dependence is due to the fact that positive a values are associated with geometrical orbits that are 
reflected an even number of times near the corner, as an optical boundary is approached, and negative 
a values with orbits with an odd number of bounces. This is explained in more detail in Appendix 
D. In the limiting case that the diffractive orbit becomes a periodic orbit, its contribution comes only 
from one of the values of a (the other cancels), and its stability matrix is M or —M depending on 
whether the number of bounces of the orbit in the corner is even or odd (cf. the discussion in 4.2). 
Thus is identical to the Maslov index of the periodic orbit in these limiting cases. 
In terms of the Presnel integral K formula (^) can written in a slightly shorter form 



Equation is the main result of this paper. It gives a uniform approximation for the contribution 
of an isolated diffractive orbit with a single corner point to the trace formula. For completeness we 
recall several definitions: (p^ = 9i — a92, where 6i and 02 are the incoming and outgoing angles at the 
diffractive point (measured from the same edge, with 9i and 02 G [0,7]) and 0^,^ is defined by 

a„„ = V2 cos(— — rirj^vl) with rig-^ = nint [ '^°" ^ — ] g ^ . (55) 
2 ' ' 27 

Note finally that the modified Fresnel function of imaginary argument (encountered when = 1) 
can be computed numerically from (|A3|). 



4 Discussion of the result 

In this section we discuss properties and the range of validity of formula ([5^). As mentioned above, 
the derivation has been done for a specific case (a wedge connected to a smooth boundary), but it 
is more generally valid because it relies only on local properties of the system near the considered 
diffractive orbit (as usual in semiclassical approximations). Hence it applies to billiards of any shape 
provided the singularity of the boundary corresponds locally to the intersection of two straight lines. 
However, the present approach has to be refined if applied to curved edges, we discuss this point in 
Appendix D. We also remind that formula ( |54|) is only valid for single diffraction. The same formalism 
can in principle also be applied to diffractive orbits with more than one diffractive point, the formulas 
become however increasingly more complex. For example, in the case of double diffraction, one has 
already 16 instead of 4 terms, and they involve also double Fresnel integrals as can be inferred from the 



treatment of diffraction at two wedges in |21|. The formulas can only be simplified if the diffraction 
in some of the corners can be treated in the GTD approximation. 

Note also, that the factor | TrM — 2o"|^^/^ in ( [5^ ) diverges for a parabolic diffractive orbit (i.e. 
when TrM = ±2) and the present approach cannot be used in this case. This is very similar to 
divergences in Gutzwiller's trace formula due to non-isolated orbits. For diffractive orbits, the fact 
that Tr M = ±2 can have several reasons, for example the diffractive orbits can appear in families as 
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is the case in a circular sector, or bifurcations of diffractive orbits can occur, or the diffractive orbit 
can become a part of a family of periodic orbits when the optical boundary is approached. The latter 
case can occur for example in triangular billiards. In this case it is however often possible to treat the 
divergent part (one of the a- values) in the GTD approximation if the diffractive orbit is well separated 
from the torus of periodic orbits, and apply the uniform approximation only to the non-divergent part 
as will be demonstrated in a numerical example in section 5. 



4.1 The GTD limit 

After these basic remarks we now study three simple limits of Eq. ( |54[ ) . The first one is the geometrical 
theory of diffraction which is valid sufficiently far away from the optical boundary. In this limit the 
argument of the K function is large and the function can be replaced by its leading asymptotic term 
in (|A4|). This immediately yields 

d^k) « - cos(fcL - z.vr/2 - 3^/4) , (56) 

vr ^8vrA;|Mi2| 

which agrees with the general formula (^) in the case of one diffractive point. Analogous formulae have 
been derived and tested in |Q, |6|, ^ . They have the advantage of allowing to treat general diffractive 



problems (other than wedge diffraction) and can easily be generalized to multiple diffraction (see (22)). 
However they diverge on the optical boundary and (as shown in the examples 5.1 and 5.2 below) they 
are unable to describe the limit that a diffractive orbit is close to become a real trajectory. 



4.2 The hmit 7 = vr/p 

Let us now study the limit that the diffraction angle 7 goes to ir/p {p G IN*). For these values of 7 
there is no diffraction since the corner can be treated by the method of images. As a consequence, 
the contributions of most diffractive orbits disappear, but there are also diffractive orbits which are 
replaced by periodic orbits which contribute to the level density according to Gutzwiller trace formula. 
Their contribution can be obtained from the diffractive contribution (|5^ ) in the limit 7 — > vr/p. The 
situation is actually slightly more complicated, since the diffractive contribution of these orbits for 
angles 7 = vr/p + e is discontinuous at e = (it changes sign) . The reason for this is that periodic 
orbits split from the diffractive orbit as e goes through zero (for example as the billiard is deformed), 
which can be considered as a kind of bifurcation. As a consequence both diffractive and periodic orbit 
contributions are discontinuous at e = 0, but their sum remains continuous. In order to discuss this 
in more detail we have to consider the cases of odd and even p separately. 

• Case 7 = ^ + £• In the limit e = the contributions from the two r/- values cancel for cr = — 1. 
The same occurs for cj = +1, except if 62 = 9i- If this condition is fulfilled one obtains 

d^{k) sign (e) r(+) dpo{k) , (57) 

where 

Mfe) = - ^°^^''^"^^^^"f^ (when 7 = -). (58) 

^ vr |TrM-2|i/2 ^ ' 2p^ ^ ^ 
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The discontinuity in Eq. (p7\j at e = is directly related to the appearance of new periodic orbits. 
This can be seen from the discussion in section 3.2: in Eq. ( |49| ) one has contributions of periodic orbits 
in the vicinity of the diffractive orbit if tta^nTa > 0, and the periodic orbits coincide with the diffractive 
orbit when tta^rj = 0. For the considered case the above inequality is equivalent to — sign(e)r+ > 0. 
Hence when e goes through zero, two periodic orbits appear (or disappear), one for each value of rj, 
assuring the continuity of the sum of contributions at e = 0. 

• Case 7 = vr/(2p + 1) + e. Now the two contributions to ex = +1 cancel as e ^ 0, and for a = — 1 
there is only a contribution if ^2 = 7 — ^i- This contribution is of the form 

d^{k) sign (e) r(_) dpo{k) , (59) 
and the periodic orbit contribution now is given by 

Mfe) = - '°'^''^"^^"^jf^ (when 7 = ^). (60) 

^ vr |TrM + 2|i/2 V ' 2p + V ^ ^ 



Comparing with (|58|), the reason for the change of sign of M is the odd number of classical 
reflections on the vertex in the case 7 = iT/{2p + 1). Generally, the stability matrix M of the closed 
trajectory in (p3) becomes equal to plus (resp. minus) the monodromy matrix of the periodic orbit 



when 7 is 7r/(2p) (resp. 7r/(2p + 1)). The explanation of the discontinuity of Eq. (59) is the same as 
above with the only difference that the condition for the existence of neighbouring periodic orbits can 
now be expressed by — sign(e)r_ > 0. 

In billiards with corners one has therefore a new kind of bifurcation: the continuity of wave me- 
chanics (in the semiclassical approximation) is not enforced by complex trajectories but by diffractive 
orbits. This effect will be demonstrated in the examples below. 

4.3 In the vicinity of an optical boundary 

The case that a diffractive orbit lies on an optical boundary, or crosses an optical boundary when the 
billiard is deformed, is very similar to the case 7 vr/p. Again the diffractive orbit contributes on the 
optical boundary at the same order of A; as a periodic orbit, but now only with half the amplitude of 
a periodic orbit. The diffractive contribution is again discontinuous since it changes sign as an optical 
boundary is crossed, and the reason for this is that a new periodic orbit arises which bifurcates from 
the diffractive orbit. More specifically, let us consider the case that for a given value of a and r] one 
has (po- — 2no-,r)7 + ijir = e where e is small. In the limit e ^ the contribution from these values of a 
and T] to the spectral density is given by 

1 1 f,\ L cosikL — ^(jTT 12) 
-r]T^ sign {e)dpo(k) where dpo{k) = , =- ■ (61) 

^ ^ ■/iTrM-2o-| 



and one can verify that the discontinuity of Eq. (61) is due to a neighbouring periodic orbit which 



coincides with the diffractive orbit at e = 0. As above, the condition for the existence of the periodic 
orbit is ao^nTa > which now is equivalent to rjeTa > 0. 
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= 2 h sin 7 



Lg = 2 h 



7 < 7t/2 




L3 = 4 h sin (7/2) L\ = = R 

7 < 7T 7 > 7T 

Figure 6: Shortest periodic and difFractive orbits in the "rounded triangular billiard" studied in Sec. 
5. For the diffractive orbit, the diffractive point is marked with a black point. The upper plot defines 
the geometry and the notations. 
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5 Some examples 

In this section we illustrate the results of the last sections with several examples. We study mainly 
a billiard consisting of a wedge of opening angle 7 whose two edges are connected by an arc of constant 
radius of curvature R. The angles between arc and wedge are chosen to be vr/S on both sides. If h 
denotes the "height" of this billiard (see Fig. 6) then R = /isin(7/2)(sin(7/2) — 1/2)""'^. This billiard 
has only one diffractive corner (at point S of Fig. 6) and the curvature ensures that the shortest 
diffractive and periodic orbits have TrM 7^ ±2 (they are displayed in Fig. 6). In the following we call 
this billiard a "rounded triangle (vr/S, 7r/3, 7)". 

For numerical convenience we restrict ourselves to angles of the form 7 = pir/q with {p, q) E IN^. 
The quantum energies are determined by expanding the wave functions around point S in "partial 
waves" which are Bessel functions times a sinusoidal function of the angle: 

nmax 

i;{r,0)= J2 Jn,/pikr)smi^e). (62) 

n=l P 

Eq. ( |62|) automatically fulfills the Dirichlet condition on the straight faces of the billiard. The 
boundary condition on the arc opposite to S is enforced in a manner identical to the improved point 



matching method presented in |32|. This results in a secular equation whose solutions are the eigen 
levels of the system. We have tested the numerical stability of our procedure by varying the number 
i^max of partial waves included in the expansion (typically Umax ~ nini [pkh/q]). For each of the 
values of 7 studied below we have computed the first 2000 eigenlevels and we have checked that they 
were determined with an accuracy of the order of 1/1000 of the mean level spacing. 



In order to visualize the importance of periodic and diffractive orbits we study in the following the 
regularized Fourier transform of the level density: 

F{x) = Vfce^^^ - «^'d(A:) dk . (63) 

Jo 

kmax is the last eigenvalue computed numerically and we take here a = 10 /k'^^^. The multiplicative 
factor ^/k in ( |63|) is included in order to cancel the singularity at /c = of a contribution of type (|5^). 
F{x) is denoted Fqm{x) if we use in ( |63| ) the exact quantum spectrum. It is denoted Fua{x) (resp. 
Fgtd{x)) when Eq. ( |5^ ) (resp. (|56|)) is used together with the periodic orbit contributions (|B21j) . 

5.1 7 near 7r/2 

This case is relevant to the discussion 4.2 above. For 7 < 7r/2 the shortest orbit is a periodic one and 
has length Li = 2/isin7. It disappears as soon as 7 > 7r/2. For 7 7^ 7r/2 one also has a diffractive 
orbit of length L2 = 2h (see Fig. 6). When 7 = 7r/2 these orbits coalesce and give a single periodic 
orbit of length 2h. Their contribution to the level density is continuous at 7 = 7r/2 as explained above: 
the contribution of L2 is discontinuous (cf. Eq. (|57|)) and this exactly cancels the discontinuity due to 
the disappearance of the orbit Li and its time reverse. 

We determined the spectrum numerically for 7 = Ttt/IB and 87r/15. The corresponding moduli 
of the Fourier transform |-F(a;)| are plotted in Fig. 7 and 8. In this figures the solid lines correspond 
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Figure 7: Modulus of the Fourier transform of the level density - see Eq. ( [6^ ) - for the rounded 

triangle (vr/S, 7r/3,7 = 77r/15). The solid line corresponds to FqMix), the long dashed line to FuAix) 
and the short dashed line to Fgtd{x) (see the text). The arrows mark the lengths of the diffractive 
and periodic orbits. The scale of lengths and wave-vectors is fixed by taking h = 1. 
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Figure 8: Same as Fig. 7 for 7 = Stt/IS. 
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to |Fqm(x)|, the long dashed lines to and the short dashed line to \FcTDix)\. The lengths of 

the included orbits are marked with back arrows. One notices the failure of GTD and the excellent 
agreement of approximation (^) with the exact result (the agreement remains equally good when 
plotting the real and imaginary part of F{x)). As stated above, it can be seen that in the vicinity of 
an optical boundary diffractive and periodic orbits contribute in the same order to the level density. 

5.2 7 near tt 

This case also pertains to the discussion of Sec. 4.2, but now in the vicinity of 7 = tt/{2p + 1). Again 
the diffractive orbit ensures continuity of semiclassical mechanics when 7 = vr: the two periodic orbits 
-L4 and -L4 of Fig. 6 disappear as soon as 7 < vr, and the contribution of L2 is discontinuous at 7 = vr, 
but the joint contribution is continuous. Here we computed numerically the levels for 7 = 197r/20 
and present the results for |-F(x)| and ?R.{F(x)} in Fig. 9. Again one can verify the failure of the 
geometrical theory of diffraction and the excellent agreement between (^) and the quantum result. 

5.3 The triangle (7r/4,7r/6,77r/12) 

In this subsection we depart from the previous examples and study the spectrum of a straight triangular 
billiard with angles (7r/4, vr/G, 77r/12) which has one diffractive wedge 7 = 7-7r/12. This billiard is of 
interest because (i) it allows to compare the performances of the uniform approximation with GTD in 
a regime where this last approximation is not inaccurate and also (ii) because it provides an example 
where our approach is not completely justified. Indeed in case of a polygonal billiard all the trajectories 
have a monodromy matrix with TrM = ±2 and this leads to a divergence in (js^). As noted in Sec. 4 
this is linked to the possible deformation of any diffractive orbit of the system considered towards a 
family of periodic orbits. Fortunately, in the present case the first diffractive orbits are far from any 
allowed family of periodic orbits and we can evaluate the if-function relevant to the divergent term 
in ( |5^ ) with the asymptotic expansion (^^): this cancels the divergence. This was done on Fig. 10 for 
the 3 first diffractive orbits of the system. 

Again the agreement with the numerical result is excellent, but here the geometrical theory of 
diffraction already gives a sensible description. Note however that the small peak due to the diffractive 
boundary orbit of length L2 is "missed" by GTD because its diffractive coefficient (p^) is zero (see the 
discussion in 2.2). The correct description of the peak was obtained by using half the contribution 
( |5^ of a usual diffractive orbit. 

For a more detailed comparison we plot the moduli of the differences \FirA{x) — FQMix)\ and 
\FaTD{x) — Fqm{x)\ in figure 11: even quite far from any optical boundary Eq. (U) supersedes the 
GTD result (|56[). This plot emphasizes the accuracy of Eq. ( [5^ ) in cases slightly out of its original 
range of application. 

6 Conclusion 

In this paper we have studied the inclusion of diffractive orbits in semiclassical trace formulae for 
billiards in which the boundary has wedge-like singularities. In many cases the simple geometrical 
theory of diffraction ^ is inadequate, especially if the energy is not very high. A consideration 
of the mathematical structure of the exact Green function near a wedge permits to remedy this 
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Figure 9: Same as Fig. 7 for 7 = 197r/20. We plot here also ^{F{x)} for illustrating the quality of 
the agreement between the phases of Fqm and F^a- 
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Figure 10: Same as Fig. 7 for the flat triangle (7r/4, 7r/6, 77r/12). The upper part displays the shortest 
orbits of the system (all three are diffractive, the classical periodic orbits occur at greater lengths). 
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Figure 11: \Fua{x) — Fqm{x)\ (long dashed line) and \Fgtd{x) — Fqm{x)\ (short dashed line) for the 
triangle (7r/4, vr/G, 7-7r/12). We consider only the three shortest orbits of the system. The following 
orbits are not taken into account and this is the reason for the increasing errors in vicinity of x « 4. 



30 



shortcoming: it leads to a uniform approximation of the Green function [16| which, in turn, aUows 
to derive contributions to the trace formula which properly account for the role of isolated diffractive 
orbits in the quantum spectrum (^). The formula was illustrated in several examples and was shown 
to give excellent agreement with numerical data. Its main feature is that it interpolates between 
the usual Gutzwiller trace formula Q and previous approaches relying on the geometrical theory of 
diffraction ggg. 

Note also - as a by-product of the present approach - that we derived a semiclassical composition 
law for Green functions for billiard systems ( Bl|jB13 ) reminiscent to the semigroup composition pro- 
perty of the propagator (see the discussion in Q). This allows to recover Gutzwiller 's trace formula in 
a simple fashion (cf. Appendix B ). Similar laws can also be obtained for the composition of diffractive 
and geometrical Green functions. 

The present work suggests further developments: (i) the result ( p4|) might be extended to allow the 
treatment of diffractive orbits in the vicinity of a family of periodic orbits; (ii) although the inclusion 
of general multiple diffraction in a uniform formula seems to be a difficult task, one may reasonably 
hope to include double diffraction in the formalism (cf. [^). (iii) Further possible extensions concern 
the treatment of other types of diffraction, like regions near curved wedges where surface diffraction 
becomes important so that contributions from creeping and whispering gallery orbits have to be 
included, or diffraction effects arising from discontinuities of the curvature of the boundary like in the 
stadium billiard. 

Finally we would like to emphasize the important role of diffraction in semiclassical approaches. 
Diffractive and periodic orbits are fundamentally different in the sense that the former are not obtained 
via a systematic h expansion in the vicinity of classical solutions of Hamilton's equations (they are 
rather linked to discontinuities of the Hamiltonian flow). However diffractive orbits provide the first 
correction to the leading order in the trace formula, with contributions typically of order smaller 
than the contributions of isolated periodic orbits. Besides, in the vicinity of optical boundaries the two 
types of orbit contribute with approximately the same order to the trace formula. An image emerging 
from our study (cf. Sec. 4.2 and 4.3) is that diffractive orbits allow to enforce semiclassically the 
continuity of wave mechanics in the vicinity of discontinuities or bifurcations of classical mechanics. 
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Appendix A The modified Fresnel function 

In this Appendix we define the modified Fresnel function K{z) used in the main text and list several 
of its properties. 

• The function K(z) (z g(D) is defined by: 
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K{z) = ^ / dy = erfc (e-^^/^z) , (Al) 

Vvr Jz 2 



where erfc is the complementary error function (see e.g. [^]). In (Al) the path of integration is 
subject to the restriction arg [y] — > a with < a < 7r/2 as y ^ co along the path, a = and 7r/2 are 
permissible if Re (iy^) remains bounded to the right. 

The function K has the following properties: K{+(X)) = 0, K{Q) = 1/2, 

K{z) + K{-z) = e-^^' , (A2) 

and 

K(f) + K{-iz) = e^^^ , (A3) 
where the bar denotes complex conjugation. 

• By successive integrations by parts one obtains the following asymptotic expansion: 



K{z) = -j= X! ( ~ ) [~^] ^^"^ \A ^ +00 and - 7r/4 < arg (z) < 37r/4 , (A4) 



2Z\fK n \2/ n 
* n=0 ^ ' " 



z 



where (l/2)„ = r(n + l/2)/r(l/2) = 1 x 3 x ... x (2n- l)/2". In the region arg {z) E ]37r/4, 77r/4[ one 
obtains an asymptotic expansion by combining Eqs. ( [A2| ) and ([A4|). 

• The interest in the modified Fresnel function comes from the following integral relation: 

hoo p— / . , \ 

dt- = 2ir7ri^ re-W4y;5^ ^ (A5) 

X t - z \ J 

where /? G IR"*", z g(D and r = sign (9 (z)). 

Hence the function ilT allows to generalize the steepest descent method to cases where poles appear 
in the integrand. As explained in the text (Sec. 2) this corresponds - in the Sommerfeld solution of 
the diffraction problem - to the occurrence of diffractive orbits near classical trajectories. We will not 
prove Eq. ( |A5|) here, it can be done easily by noting that (t — z)~^ = it /q^°° exp[zT(z — t)x\dx (cf. the 
evaluation of integral (C6) in Appendix C|). 



Appendix B Composition law for Green functions 

In this Appendix we derive a simple semiclassical composition law for Green functions which is ex- 
pressed by integrals over the boundary dB of the billiard. Although the formulae established below 
are simple and natural from the point of view of Balian and Bloch's multiple reflexion expansion 
to our knowledge they have not been clearly stated in the literature. The composition law can be used 
in order to simplify expressions obtained from the boundary element method (cf . (p|) ) . We first prove 
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the semiclassical version and then give the exact formulation of this law due to Balian and Bloch. We 
further show that it allows to derive Gutzwiller's trace formula in a straightforward manner. 

We assume in this Appendix that the boundary dB is smooth everywhere. The semiclassical 
version of the composition law has the form 



(-2)" / dsi...dsnGoin,r\E)dn,Goir2,n,E)...dn^Goif,fn,E)^Gt\r,f',E) . (Bl) 

where rj = r(sj). The approximate sign signifies that the evaluation is done by approximating the free 
Green function Go by its leading asymptotic term for large argument and evaluating the integrals in 
stationary phase approximation. The function G^\f, r E) on the right-hand side of Eq. ( piD is the 
part of the semiclassical Green function from all trajectories with n bounces on the boundary between 
r ' and f 

Gt\r,r',E) = ^ \ expjiMW - - i^} . (B2) 

Here ^'^"^ denotes the length of the trajectory, the number of conjugate points from r ' to r 

plus twice the number of reflections on the boundary, and M^") is the stability matrix for unit energy. 
An index ^„ of the above quantities has been omitted in order to simplify the notations. 

Eq. is proven by mathematical induction. For n = it is correct since 

Go (f , f', E) « G(o) (f , f, E) , (B3) 

and one has to show that 

/ := (-2) / dsiGfXri,r\E)dn, Go{r,n,E) « G'i::+'Hf,r' , E) . (B4) 

We will use the following notation at a point fi of the boundary: primed quantities correspond 
to the incoming trajectory and unprimed quantities to the outgoing trajectory. We will denote the 
momentum of a classical trajectory by a vector p of modulus k whose direction is the direction of 
propagation of the classical particle. The momentum of an outgoing trajectory is pi and Oj is the 
angle between the normal vector hi of the boundary (which points outside) and —pi. The momentum 
of an incoming trajectory is p^, and a[ is the angle between fii and p^. For this choice a'^ and a, both 
lie in the interval between —it/2 and 7r/2. In terms of the local coordinate systems of the trajectories 
with coordinates parallel and perpendicular to the trajectory, the tangential and normal vectors of 
the boundary can be written as: 

hi = — cos Oi eii + sin Oi e± = cos a- eji — sin a'i e\ , 



ti = — sinoj eii — cosoj e_L = sin a • Cn + cos a' e 



We continue by evaluating the integral in Eq. ( |B4| ) using the stationary phase approximation. The 
normal derivative of the Green function is given in leading semiclassical order by 

dfi. Go (fi+i ,fi,E) K. -ihi- Pi G^J (f^+i , fi,E) = ik cos ai G^^ (fj+i ,fi,E) , (B6) 
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and in ( [B4D the stationary points are determined by the condition 

= A[;(0)(^^^^) + zW(f,,f')] = £• [-^ + ^] = sinai + sina^ , (B7) 
asi k k 

i.e. by ai = —a'l, which is the condition for elastic reflection. The sum over all stationary points 
thus expresses the integral / by a sum over all trajectories with n + 1 reflections on the boundary. In 
Eq. ( p37| ) and in the following the length is given two arguments when it is necessary to specify the 
starting and end point of the trajectory. 

For the determination of the second derivatives of the lengths at a boundary point Si one has to 
evaluate the derivatives of the angles ai and a- which consist of two parts. One from the change of 
the normal vector with Sj and one from the change of the direction of the trajectories. 



d^l^''\run-i) _^^^,da[ 1 . cos a[ dp'i\ _ cosa^ cos^ M 
cos a,- -J— — cos a,- — ^ i r tt— — d 1 



d^l^^-^fi+ufj) _„^^^dai _„„^^Y 1 cos aj dpi A _ cosa^ , cos^ M^^ 



-,2 „/ nAn) 



H ^"22 



'12 



(B8) 



At a stationary point it follows with these relations that 

where M("+^) = M(°) Bi M(") and the matrices M(°) and Bi are given by 

M'»' - ( ' ). B,J M . (BIO) 



\ / \ if i COS ai 

The matrices M and -B correspond to the linearized flow near the considered trajectory. Note 
that our definition is slightly different from usual conventions (see e.g. [Q, considering that here 
IpI = k, the Mi2 (resp. M21) matrix element would be generally divided (resp. multiplied) by k. Here 
we work with the stability matrix at unit energy: this choice is connected to the scaling property of 
the dynamics in billiard systems. It does not affect the trace and the determinant of the matrix and 
allows to have energy independent matrix elements with a simple geometrical meaning. 

Now the stationary phase approximation for the integral in Eq. (|B4|) is carried out and results in: 



^ cosai exp{iA;Z('^+i) -i^i/^")} /• , , , cos^ ai m}" 
/ , =^ / dsi exp{— zfc TTTT — 



V ^ ^ exp{ifc/("+i) - z-z.("+i) - i—} , (Bll) 

where = + and 
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^ ^(n) + 2 + I ^ sign(M;r'J) = sign(MgJ) , 

\ if sign(M}^+^)) / sign(Mf^)) . ^ ' 

Eq. ( |B12| ) coincides with the expected definition of the Maslov index: the number of 

conjugate points from r ' to r plus twice the number of reflections on the boundary; an additional 
conjugate point has occurred between fi and r iff sign(M}2'''^^) = sign(M^2^) (remember that there 
is one sign change due to the reflection on the boundary). This completes the proof of Eq. ( [B4D and 



thus also of Eq. (Bl). 

A further relation follows from the fact that the evaluation of the integral in Eq. ( |BlD does not 
depend on the order in which the stationary phase approximations are carried out. Thus one can 
conclude directly that 

(-2) I ds^G^rXri,r\E)dn,G^:^\T,n,E) ^ Gt^^+^\r,r\E) . (B13) 
JdB 

Equations and (^) were derived in the semiclassical approximation by evaluating the boun- 



dary integrals only locally in the vicinity of stationary points. For that reason the same composition 
law can be applied in order to obtain the contributions of the geometrical orbits in billiards with 
corners; this is done in Eq. pO|). 

We note that equations (|B1|) and ( [B13| ) are the semiclassical versions of exact relations for the 
Green function G of a billiard system. These exact relations are obtained by a multiple reflection 
expansion of the Green function G [^: 

oo 

G(f, f E) = Y, (f, f E) , (B14) 

n=0 

where 

G(")(f,r",£;) = (-2)" / dsi . . .dsnGo{ri,f' ,E)dn,Ga{r2.ri,E) . . .dn^Goir^rn.E) , (B15) 
JdB 

and the equation analogous to (B13) follows directly. 



Finally, we show that Gutzwiller's trace formula can be obtained in a straightforward way by using 
(pil). From the boundary element method one obtains 

-1 oo -| 7 

d{k) = d{k) + -9 V -— TrQ"(A:) , (B16) 
vr , n dk 

n=l 



where 



TrQ"(fc) = (-2)'^ / dsi . . .dsndn,Go{r2,ri,E)dn,Go{rz.r2,E) . . .dn^Go{ri,rn.E) . (B17) 
JdB 



I OB 

With Eq. (pl|) it follows that 
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IT dk , n Job 

n=l 



d(k)--^y-^^^l^^ I d.exp{iM(«"i)(f,f)-i-i.("-i)-i;}. (B18) 



The stationary phase condition is again given by sin a = — sin a' and thus the integral yields 
contributions from periodic orbits with n specular reflections on the boundary. More accurately, it 
gives n/rpo (identical) contributions for every periodic orbit where Vpo is the repetition number of the 
orbit, since there are n/rpo different starting positions r = r' on dB. 

The derivatives of the angles a and a' now have additional contributions since both initial and 
final points of the trajectory are changed by varying s 
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ds 
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M(n-l) . / 
II + COS a 



(B19) 



M 



(n-l) 
12 



It then follows at a stationary point that 



ds'- 



cos^ a (Tr 



po 



po jl2 



(B20) 



where = BiM^''-^\ 



and the stationary phase approximation results in 



d{k) = m + - E E 



/{") 

''po 



n=i in,po rpo Y |Tr Mj 



^po 



:COs{fc4") - 



"po 



2 t^po i 1 



(B21) 



where 



H-po 



+ 2 + 



if 
if 



^po 

po , 



i2/(Triv4^) 



2) > 
2) < 



(B22) 



Note that the derivation presented here has the same starting point as Ref. |^^. But the composi- 
tion law (|B1| ) permits to bypass the computation of large determinants of |2^. Furthermore, it allows 
to keep track of the Maslov indices (which were not derived in |2^ ) in a simple way. 

Finally we add a remark on ghost contributions. In general, the semiclassical approximation for 
the Green functions G^"'\f,'r',E) can also contain contributions from ghost trajectories that satisfy 
the stationary phase conditions but have parts that are outside the billiard region. These ghost orbits, 
however, do not give a contribution to the level density d{k) since they cancel with ghost contributions 
from different n or from d{k) 28, 2£, 34 1. 
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Appendix C Evaluation of a diffraction integral 



In this appendix the integral 

■ 2 • 2 

/oo rioo filCiS — tCZ 
ds / dz , (CI) 
-oo J —ioo Z -\- S Sq 

is evaluated for positive c and real non- vanishing a and sq. This is the basic integral which appears 
in the derivation of the uniform approximation for diffractive contributions to the trace formula. 

First the z-integral is rotated onto the real axis. The rotation is performed counter-clockwise: 
since c > this yields no contribution from infinity. There are, however, poles of the integrand on 
the real z line. We take them into account by giving to sq a small imaginary part sq — > so + ^^o^ and 
consider the limit e — > in the end. Here gq = sign(so) and e > 0. For this choice one obtains a pole 
contribution from the rotation of the z-integral for those values of s for which (sq — s) has a different 
sign than sq. One obtains 

■ 2 ■ 2 

p — ICZ oo p — ICZ . > ^2 

dz — = lim / dz + 27ricJo e"'^^^ " G(ao(s - sq)) • (C2) 

-ioo 2; + S - So ^-^0 Joo Z + S - So - lOnE 



We consider now the two contributions of the r.h.s. of Eq. (C2) to the integral in (CI) separately, 



/ = /o+/i, where /q contains the pole contribution and l\ the contribution from the rotated z-integral. 
For 7o we have 



/oo 
ds G(o"o(s — So)) 6 
-oo 



ias^ — ic{s — so)^ 



27riao /"cise^«^'-^^(^-|^ol 



, exp{-::^} erfc{ 7'^°' } , (C3) 

\j—i[a — c) a — c y'—i[a — c) 

where erfc is the complementary error function (see e. g. [^]). Ii has the form 

• 2 ■ 2 

/oo POO pZClS — ICZ 
ds / dz . (C4) 
-oo J-oo z + s - So - iaos 

By a linear transformation of the variables 

cz as 

u = z + s 1 V = 1 , (Co) 

a—c a—c 

the double integral splits into a product of two single integrals 

. ac 9 



hm / d7;e^(«-^)^ / du^—— . (C6) 

e^Oj_oo J-oo U-SQ-iaQE 
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In ( |C6| ) the integral over v can be computed easily. Furthermore, the denominator in the u-integral 
can be expressed in terms of an integral 



lim 

-i[a — C) e^O 



du 



ac 2 
-I u 



dwe a - c 



iaow{u — So — idQe) 



IT 



' vr a 



■i{a 



lac 



dw e 4ac 



Co -11 

W + 1\Sq\W 



y/—i{a — c) ' ~ a — c 
The whole result I = Iq + Ii is given by 



exp{ — ^;£fO| erfc{— i|so| 



lac 
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exp{- 



tacs, 



} 



erfc{ 



-ia\S()\ 



^J-i{a - c) 



} — erfc{— z|so| 



lac 



(C7) 



(C8) 



It is convenient to rewrite this result in a form in which the phases of the complex arguments 
of the error functions are always between — 7r/2 and 7r/2. This can be done by considering all the 
possible cases for the signs of a and (a — c) separately and using the relation erfc(2:) = 2 — erfc(— z). 
The results for the different cases can be combined again and written in the form 



vr 



where Ua = sign(a) and r 



sign(a/ (a — c)). 



erfc{ 



|aso| 



:} - erfc{|so| 



ac 



(C9) 



Appendix D Curved wedges 

In this section we discuss the effect of curved wedges on the contributions of diffractive orbits to the 
level density. The uniform approximation ( |5^ ) has been derived for a boundary with zero curvature 
on both sides of the corner. It has to be modified for curved wedges, otherwise the sum of diffractive 
and periodic orbit contributions is not continuous any more as an optical boundary is crossed. Addi- 
tional complications can arise due to surface diffraction effects, i. e. creeping orbit or whispering orbit 
contributions can interfere with the diffractive orbit contributions. We will discuss when these effects 
have to be taken into account, but we will modify the uniform approximation only in those regions in 
which these additional effects can be neglected. 

The modified formula is derived by using a method of Ref. [^] for obtaining a uniform approxima- 
tion for the Green function in the case of a curved wedge (see also Ref. |jl^ ) . We refer to the original 
references for a discussion of this method and state here only the result which consists of a change 
of the argument of the Fresnel function in (21) such that the approximation is continuous across an 
optical boundary. For the diffractive orbit contribution to the level density, this has the consequence 
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53,54): there are additional contributions to 



that only the stability matrix M is changed in Eqs. ([ 
M from reflections on the curved boundary. 

In order to discuss these modifications, we first list several properties of geometrical orbits cor- 
responding to an optical boundary which is specified by the values of a, rj and ria^r)- In particular, 
we consider the trajectories which contribute to the Green function and list for them the number of 
reflections on the boundary and the side of the corner on which the first reflection occurs. Furthermore, 
we give restrictions for the numbers rio-^r; which are implied by their deflnition. 



a = +1, rj = +1: rifj^n > 0, {2na,n) reflections, first on the line 6 
+1, r] = —1: n(j,r; < 0, (— 2no-,r;) reflections, first on the line 



7- 



a 



0. 



• a = —1, 7] = +1: n^^r] ^ 0, (2no-,r) — 1) refiections, first on the line 9 = 'j. If rio-^r; = the optical 
boundary cannot be reached. This case can occur only for 7 > vr. 

• a = —1, T] = —1: n^^j^ < 1, (1 — "^na^r]) refiections, first on the line 9 = 0. The optical boundary 
cannot be reached if n^^^j^ = 1. This case can occur only for 7 > vr. 

With these properties we can now discuss the modification of the stability matrix M in the case of 
curved wedges: M then has an additional contribution for every of the refiections mentioned above. 
In the following we denote the limits of the radii of curvature as the corner is approached from either 
side by Rq and where the first one corresponds to the side ^ = and the second one to 9 = j. 
Then M has to be replaced by BM where 



and 



E 



E 



j=0 

E 
E 

j=0 



Rosm{2n-9i] 

2 

Rosm{2n + 9i] 

2 

Rosm{2n-9i] 

2 

Rosm{2j-f + 9i] 



B 



na-,ri — l 

+ E 

i=o 

+ E 
+ E 

j=0 

+ E 




(Dl) 



i2^sin((2j + 1)7-^1) 

2 

i?^sin((2j- 1)7 + ^1) 

2 

i?^sin((2i + 1)7-^1) 

2 

R^sm{{2j -l)j + 9i) 



if a = +l,r/ = +1 
if cr = +l,r/ = -1 
if a = -l,r/ = +1 

if 0" = — 1,7/ = — 1 



(D2) 



This approximation is only valid as long as all sine-functions in (^^) are positive and not close to 
zero. The case of an almost vanishing sine-function corresponds to almost grazing incidence on a side 
of the corner. Then surface diffraction effects become important and interfere with the diffractive orbit 
contribution, and the uniform approximation is no longer valid. In the case that some sine-functions 
are negative and not small, the orbit is not close to an optical boundary and the GTD-approximation 
can be used (it is the same as in the case of non-curved wedges). 
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There is a disadvantage of the definition of B given above. Since there are two possibihties for 
choosing 9i and 02 corresponding to the two arms of a diffractive orbit in a corner, it follows from ( p2| ) 
that the uniform approximation is not uniquely defined. (Note that both cases have to be checked for 
deciding whether surface diffraction effects are important.) 

The non-uniqueness of the approximation is a direct consequence of the fact that the uniform 
approximation for the Green function of Ref. [p!^ ] is not symmetric under 9i ^ 62. It is another 
example for the non-uniqueness of uniform approximations (cf. the discussion in Sec. 2.3). However, 
as an optical boundary is approached both choices give the same result as they should. Let us discuss 
in more detail the difference between these two choices. It can be shown that interchanging 9i and 
02 amounts to evaluate (D2) with 0i = a02 + 2n„^rjl — 'H'^ instead of 0i. This then directly suggests a 



possible way by which this ambiguity can be removed, namely by replacing 0i in ( |D2| ) by the average of 
both values which is {0i+a02 + 2nf^^fi'y — riTT)/2. As an optical boundary is approached this combination 
again becomes identical to 0i. 
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